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ABSTRACT 

This  paper  is  a  continuation  of  [ 2] .  Here  we  apply  the  ideas  of  that 
paper  to  the  linearized  shallow  water  equations ,  or,  equivalently,  the 
linearized  two-dimensional  isentropic  Euler  equations  of  gas  dynamics. 
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SIGNIFICANCE  AND  EXPLANATION 


When  one  computes  a  numerical  approximation  to  the  solution  of  a  partial 
differential  equation,  it  is  sometimes  necessary  to  restrict  the  computation 
to  a  domain  which  is  only  a  portion  of  the  domain  on  which  the  problem  is 
naturally  defined.  This  is  done  when  one  is  interested  in  the  behavior  of  the 
solution  only  on  that  part  and  when  a  computation  over  the  entire  domain  would 
be  prohibitively  expensive. 

In  such  cases  a  portion  of  the  boundary  of  the  computational  domain 
represents  merely  the  edge  of  the  computation  and  corresponds  to  nothing 
physical.  The  task  of  finding  suitable  boundary  conditions  to  impose  on  this 
part  of  the  boundary  can  present  some  analytical  and  numerical  difficulties. 
This  work  is  concerned  with  one  such  difficulty  which  arises  in  limited-area 
meteorological  computations.  It  is  also  related  to  the  problem  of  finding 
conditions  which  prevent  the  nonphysical  reflection  of  waves  at  an  artificial 
boundary. 
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BOUNDARY  CONDITIONS  FOR  SUPPRESSING  RAPIDLY  MOVING 
COMPONENTS  IN  HYPERS^LIC  SYSTEMS  lit  AN  EXAMPLE 

Robert  L.  Higdon 

1.  Introduction 

This  paper  continues  the  discussion  of  hyperbolic  systems  begun  in  [  2] . 
The  goal  of  this  work  is  to  find  boundary  conditions  which  prevent  rapidly 
moving  high-frequency  waves  from  entering  the  spatial  domain  on  which  the 
system  is  defined,  in  the  present  paper  we  assemble  the  ideas  of  [2]  and 
apply  them  in  detail  to  an  example  of  physical  interest.  This  is  done  partly 
for  the  sake  of  the  example  and  partly  for  the  sake  of  developing  the  concepts 
of  [2]  into  a  definite  method. 

We  will  consider  the  system 


The  unknowns  u,  v,  p  are  functions  of  x,  y,  and  t.  For  most  of  this 
paper  the  spatial  domain  we  consider  will  be  the  half-space  x  >  0.  We  assume 
«  f  0  and  |a|,  |$ |  «  c.  The  latter  assumption  guarantees  the  presence  of 
propagation  speeds  having  substantially  different  magnitudes. 

The  system  (1.1)  is  a  symmetric  version  of  the  linearised  shallow  water 
equations.  In  this  notation  (-a,-0 )  is  the  velocity  of  the  flow  about  which 
the  system  has  been  linearized,  and  c  is  the  speed  of  gravity  waves  relative 
to  the  underlying  flow. 

Research  performed  at  Stanford  University  with  partial  support  from  the  Office 
of  Naval  Research  Contract  N00014-75-C-1132.  Written  at  the  Mathematics 
Research  Center,  University  of  Wisconsin  -  sponsored  by  the  United  States  Army 
under  Contract  No.  DAAG29-80-C-0041.  This  material  is  based  upon  work 
supported  by  the  National  Science  Foundation  under  Grant  No.  MCS-79  2706  2,  Mod.  1. 


Up  to  multiplicative  factors/  the  dependent  variables  u  and  v  are  the 
perturbations  in  the  components  of  velocity.  The  variable  p  is  the  pertur¬ 
bation  in  the  geopotential.  The  coefficients  in  the  undifferentiated 

term  are  due  partly  to  Coriolis  effects  and  partly  to  the  process  of 
linearization.  \ 

The  linearized  two-dimensional  isentropic  Euler  equations  of  gas  dynamics 
can  also  be  written  in  the  form  (1.1).  In  this  case  c  is  the  speed  of  sound 
waves,  and  p  is  essentially  the  perturbation  in  the  density. 

In  order  to  find  the  desired  boundary  conditions,  we  first  transform  the 
system  to  an  approximate  diagonal  form  in  such  a  way  that  each  of  the  new 
dependent  variables  can  be  identified  as  a  slow,  incoming  fast,  or  outgoing 
fast  portion  of  the  solution.  We  then  set  the  incoming  fast  part  equal  to 
zero  at  the  boundary,  at  least  approximately.  The  calculations  presented  here 
follow  the  outline  given  in  Section  5.4  of  [2]. 

Pseudo-differential  operators  are  used  extensively  throughout  this 
work.  The  use  of  these  operators  is  analogous  to  the  formal  manipulations  of 
Fourier  transforms  which  are  used  to  state  the  outline  in  [ 2] •  The  operators 
are  introduced  in  order  to  treat  systems  with  variable  coefficients  and  to 
settle  a  technical  point  mentioned  in  Section  3.2  of  [2].  An  informal 
discussion  of  the  operators  is  given  in  Section  4.1  of  [2],  and  in  the 
Appendix  of  the  present  paper  we  define  them  precisely  and  state  some  of  their 
basic  properties.  In  this  work  they  are  regarded  as  a  formalism  for  studying 
the  high-frequency  asymptotic  behavior  of  the  solution. 

In  Section  2  we  uncouple  the  system  (1.1),  and  in  Sections  3  and  4  we 
derive  boundary  conditions  from  the  results  of  the  uncoupling.  Numerical 
tests  of  the  conditions  are  discussed  in  Section  5.  An  analytical 
justification  of  the  boundary  conditions  is  outlined  in  Section  4.2  of  [2]. 


2.  Uncoupling  tha  gygtfi 


Ha  first  diagonalize  tha  coafficiant  tutrix  of  tha  normal  derivative 
in  ( 1 . 1 ) .  To  do  this  we  use  tha  matrix 


(2.1) 


i  /  £  1  1 

—  1/2  0  0 
✓  2\  0  1  -1 


Tha  columns  of  this  matrix  are  normalised  aiganvactors  of  tha  coafficiant 
matrix  in  question.  When  wa  multiply  (1.1)  on  tha  left  by  tha  inverse  of 
(2.1)  and  make  tha  appropriate  change  of  dependant  variable,  tha  result  is  tha 
system 


(2.2) 


where 


w  +  D*  +  Ew  , 
x  y 


W  -  I  —  (u+p) 

(2.3)  /  2 

\/2 (u  p)y 

D  and  E  can  be  written  explicitly,  but  we  will  wait  until  we  solve  for  wx. 

The  system  ( 2. 2)  can  be  used  to  make  a  preliminary  identification  of  the 
slow,  incoming  fast,  and  outgoing  fast  portions  of  the  solution.  Since 
|a|  «  c,  we  can  say  roughly  that  the  first  component  v  is  a  slow  component 
and  that  u+p  and  u-p  are  incoming  and  outgoing  fast  components, 
respectively.  In  order  to  suppress  the  incoming  fast  part  of  the  solution  one 
could  therefore  require  u+p  ■  0  at  the  boundary  x  *  0.  If  a<0,  then 
one  would  also  prescribe  a  value  for  v  in  order  to  have  a  well-posed 
problem.  The  trouble  with  this  approach  is  that  it  ignores  the  effect  of  the 
terms  Dw^  and  Ew.  Our  identification  of  the  various  parts  of  the  solution 
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In  analogy  with  the  outline  In  Section  5.4  of  [2),  we  write  the  system 


(2.4)  In  the  form 


(2.6)  wx  ■  Hw  +  Cw  +  Ew  , 
where  H  la  the  operator  whoae  symbol  la  given  by 

(2.7)  0H  -  i?A  +  lrn+B  , 

and  B  la  the  operator  with  aymbol  iw(1-$)B.  ♦  la  the  eaooth  cutoff 

function  which  la  used  to  reatrlct  attention  to  the  conical  aet  T  in  the 
(«.£)  apace  in  which  faat  waves  can  be  found.  As  noted  In  the  outline,  the 
term  Bw  can  be  regarded  as  an  insignificant  error  term. 

In  order  to  uncouple  the  leading  order  part  of  ( 2.6 ) ,  we  need  to  find  a 
symbol  q  such  that  q3„q  1  la  closer  to  diagonal  form  than  o  ,  at  least 

n  n 

4) 

for  small  £-.  To  do  this  we  will  use  the  perturbation  method  given  in  the 
Appendix  in  [2].  That  is.  we  will  find  a  matrix  M  such  that 

(I  +  |  ♦M)(A  ♦  ♦BMl  ♦  |  ♦M)“1 


(2.8) 


-  diagonal  matrix  +  (o(j|r  $)2) 

We  will  then  let  q  -  I  +  jjr  *M  and  apply  to  (2.6)  the  pseudo-differential 
operator  which  is  defined  by  the  symbol  q. 

Satisfying  the  condition  (2.8)  amounts  to  setting  the  off-diagonal 
elements  of  MA  -  AM  +  B  equal  to  zero.  A  calculation  shows  that  M  can  be 
taken  to  be 


(2.9) 


1 


0  -(o-c)  -(c+c) 


M  -  —  |  a 

/2  la 


0 

0 


-1 


and  that  (I  +  6M)(A  +  eb)(I  +  em)  is  then  equal  to 


2 

a 


(2.10) 


1 

a-r 


3 

a 


+  e 


n 


+  o(e2) 


The  off-diagonal  elements  in  (2.9)  are  determined  uniquely  by  the  condition 
(2.8)/  but  the  diagonal  elements  may  be  chosen  arbitrarily.  For  convenience 


we  have  set  these  equal  to  zero. 

Mow  define  the  symbol  q  by 
(2.11)  q  -  I  +  |  +M  / 

where  M  is  given  in  (2.9),  and  let  Q  be  the  pseudo-differential  operator 
with  symbol  q.  When  the  operator  Q  is  applied  to  the  system  (2.6)#  the 


result  is 


(2.12)  (Qw)x  -  (QHQ'bQw  +  (QCQ- 1  +  Q^^Qw  QEw  . 

Here  Qx  denotes  the  operator  with  symbol  qx/  and  Q~  ^  is  a  parametrix,  or 
approximate  inverse/  of  the  operator  Q.  Par ametr ices  are  defined  in  the 
Appendix. 

The  composition  law  for  pseudo-differential  operators  implies  that  the 
leading  symbol  of  QHQ-1  is  equal  to  qU„q  1 .  According  to  (2.7)  and  (2.11) 
this  is  given  by 


qCrHq-1  -  (I  +  |  ♦M)(i?A  +  id*B)(I  +  |  ♦M)“1  . 

(it 

If  we  factor  out  it  and  identify  p-  $  with  e  in  ( 2. 10 ) ,  we  can  conclude 

2  * 

-1  oj  .  2 

that  qo  q  is  equal  to  cr  +  (p—  ♦  ),  where 

n  G  s 


(2.13) 


°G  "  15 


1 

a-c 


+  iw$ 


1 

a+c 


Let  G  denote  the  operator  whose  symbol  is  o  . 

G 


be  written  in  the  form 


The  system  ( 2. 1 2)  can  then 


(2.14) 


3  w 

0 

3  x 


Gw„  +  Zw_  + 


u  2 
o(-p  Ow 


+  o(u>(1-$))w 


9 


where  Wq  ■  Qw. 


The  operator  Z  is  associated  with  the  zero-order  terms  appearing 
explicitly  in  ( 2. 1 2)  and  with  the  terms  of  order  zero  or  less  which  arise  in 
the  expansion  of  the  symbol  of  QHQ- 1 .  in  a  moment  we  will  discuss  this 

u2  9 

further.  The  term  o(|—  $  )w  denotes  the  effect  of  an  operator  whose  symbol 

2 

is  dominated  by  ♦  and  which  is  a  result  of  the  error  in  uncoupling  the 
leading  order  part  of  the  system  (2.6).  The  term  o  (w(1-$))w  represents  the 
effect  of  the  operator  E  which  appears  in  (2.12).  This  is  an  insignificant 
error  term. 

2.2.  Uncoupling  the  term  of  order  zero 

The  system  is  now  partially  uncoupled  near  normal  incidence/  since  the 
symbol  of  6  is  diagonal.  We  next  need  to  reduce  the  coupling  caused  by  the 
zero-order  operator  Z,  which  is  given  by 

Z  -  QCQ-1  +  QXQ~ 1  +  (terms  of  order  zero  or  less 

(2.  IS) 

arising  from  the  expansion  of  QHQ- 1 )  . 

The  immediate  task  is  to  identify  explicitly  the  leading  symbol  of  this 
operator. 

A  short  calculation  shows  that  the  first  two  terms  in  (2.15)  are  given  by 
QCQ-1  +  QXQ_1  -  C  +  order  (-1)  +  o^  $)  . 

(A)  (I) 

We  will  regard  terms  of  order  ■=■  ♦  as  error  terms ,  since  j  ♦  is  no  larger 

2  *  4 
0)  2 

than  the  term  $  which  has  already  appeared  in  (2.14).  The  terms  of 
order  zero  or  less  in  the  expansion  of  QHQ- ^  can  be  calculated  using  the 
composition  law  for  pseudo-differential  operators.  The  calculation  is  long 
and  tedious  and  will  not  be  given  here.  Details  can  be  found  in  [1].  The 

♦  M  +  o(jr  ♦)  +  order  (-1) 


terms  are  given  by 


Terms  involving  derivatives  of  $  should  also  appear  here ,  but  these  have 
been  omitted  since  they  are  nonzero  only  near  the  edge  of  T .  The  symbol  of 
the  operator  Z  is  then  given  by 

C  +  $MA  +  o(tt  ♦)  +  order  (-1)  . 

y  5 

The  system  (2.14)  can  therefore  be  written  in  the  form 


(2.16) 


3w0  U>2  u 

-  Qw0  +  ZQw0  +  o(^-  *>  +  o(|-  *)w 


+  o(w(1-$))w  +  (order  -1)w  , 


where  zQ  is  the  operator  whose  symbol  is  $(C  +  MA  ).  The  matrix  C  has 
been  split  into  the  sum  $C  +  (1-$)C  in  order  to  give  a  neater  form  to 
certain  formulas  which  will  appear  later. 

We  now  use  the  method  outlined  in  Section  5.4  of  [2]  to  uncouple  the  term 
of  order  zero.  That  is,  we  apply  an  operator  of  the  form  I  +  K  to  (2. 16), 
where  K  has  order  -1,  and  then  make  the  corresponding  change  of  dependent 
variable.  We  then  choose  K  so  that  the  zero-order  operator  in  the 
transformed  system  has  a  diagonal  leading  symbol.  A  short  calculation  shows 
that  this  operator  is  KG  -  GK  +  Zg.  The  calculation  requires  the  expansion 
(I+K)  1  ~  I-K  +...  ,  which  is  given  in  the  Appendix.  It  follows  that  we  need 
to  choose  K  so  that 


(2.17) 


o„o  -  oo  +  $ (C  +  MA  )  ■  diagonal  matrix  . 
K.  Cj  v»  k.  y 


Here  o  and  o  are  the  symbols  of  K  and  G,  respectively,  and  the  third 

K  G 

term  is  the  symbol  of  ZQ.  0Q  is  given  in  (2. 13),  A  and  C  are  given  in 
(2.5),  and  M  is  given  in  (2.9). 

A  method  for  solving  (2.17)  is  given  in  the  Appendix  in  [2].  After  a 
certain  amount  of  labor  one  can  obtain 


( 2.18) 


°K“i5d  -j*)c  l*2; 


t 
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whara 


*12  “  -  •<VeT,! 

k,3  mjz  [<-Y21^a)(a+o)  ♦  a<«y«y» 

k2’  mj=  <-a(T12*Y32>  -  V°‘C)1 
Si  *  IalT12’T321  +  ay(a+e)I 

ka  ’1  <J»1inl3^31<1,33,“BO» 


S2-1  (Sl,1rlJ^3rT33,(a-C’ 


Equation  (2.17)  doas  not  impose  any  conditions  on  tha  diagonal  elements  of 
K.  For  convaniance  we  have  sat  these  equal  to  zero. 

The  operator  I  +  K  transforma  the  system  (2.16)  into  the  form 


3w. 


■  Gw1  +  (diagonal  operator  of  order  zero)w^  +  (order  (-1))w 


(2.19) 


+  o(^—  $  )w  +  o(|-  $)w  +  o(u>(1-$))w  , 


where  w1  «  (I+K)Wq  *  (I+K)Qw.  The  symbol  of  K  is  given  in  (2.18),  the 
symbol  of  Q  is  given  in  ( 2. 1 1 ) ,  and  the  components  of  w  are  given  in 
(2.3).  This  represents  all  of  the  uncoupling  which  we  will  do  for  this 
system. 
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3.  Derivation  of  boundary  conditions  for  the  fast  part  of  the  solution. 

It  is  now  time  to  use  the  results  of  the  uncoupling  process  to  derive 
boundary  conditions  for  the  system  (1.1).  it  is  necessary  to  identify  the 
incoming  fast  component  for  the  partially  uncoupled  system  (2.19)  and  then 
find  conditions  which  suppress  this  component  at  the  boundary  x  «  0. 

The  symbol  of  the  operator  G  which  appears  in  (2.19)  is  given  in  (2.13) 
and  is  equal  to 


Since  |a±c|  >>  |a  | ,  the  second  and  third  components  of  w^  in  (2.19)  are 
the  rapidly  moving  portions  of  the  solution.  The  second  is  the  incoming 
component,  since  a-c  <0.  We  need  to  use  the  identity  *  (I+K)Qw  to  find 
an  explicit  formula  for  this  component  and  then  use  this  formula  to  find 
suitable  boundary  conditions  for  (1.1). 

The  dependent  variable  w^  is  given  by 

w 1  -  (I+K)Qw 

(3.1) 

-  "(I  +  Ofc)"  #  "(I  +  |  ♦Mj-w  . 

Here  we  have  used  the  expression  ( 2. 1 1 )  for  the  symbol  of  Q.  The  quote  marks 
denote  pseudo-differential  operators  having  the  stated  symbols,  and  the  small 
circle  denotes  composition  of  operators.  (3.1)  can  be  written 

(3.2)  v  -  -(♦  +  0k)"  0  "(♦  +  |  Ww  +  o(  1— $ )  . 

In  order  to  produce  cleaner  formulas  later  on  we  have  used  the  cut-off 
function  $  to  restrict  the  solution  to  the  set  T  in  the  (m,()  space  in 
which  the  fast  waves  can  be  found. 
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According  to  the  composition  law,  the  symbol  of  the  composition  in  (3.2) 


is 

(♦  +  ♦  jF  ♦*>  +  X  +  ®k>  <♦  ♦  f  ♦*) 

(3.3) 

*  i  h  <♦  *  V  Ip  {*  *  f  *M>  +  order  (_2)  * 

The  derivatives  of  ♦  +  with  respect  to  w  and  £  are  of  order  -2, 

since  o  has  order  -1  and  we  are  ignoring  derivatives  of  $ .  (3.3)  is 

& 

therefore 

<♦  +  aKW  +  jp  ♦**)  +  order  (-2) 

(3.4)  -  *2  +  |  +  |  ♦o^M  +  order  (-2) 

-♦2  +  ^K  +  r+2M  +  °(^2  *2)  +  order  (‘2)  * 

To  obtain  the  last  line  we  used  (2.18)  to  conclude  <JR  -  o(^-  <t> )  >  (3.2)  and 

(3.4)  now  imply 

(3.5)  Wl  -  "t*2  +  ♦ojc  +  |  *Si)"w  +  o(~  ♦  *)  +  o(1-4)  +  order  (-2)  . 

The  error  terms  in  (3.5)  can  be  ignored  in  the  system  (2.19),  since  their 


a  is  given  in  (2.18)  and  M  is  given  in  (2.9). 

K 

The  second  component  of  z  is  the  incoming  fast  component  which  we  need 
to  suppress.  According  to  (3.7)  and  the  expressions  for  0R  and  M,  this 
component  is 

*2(x, y*t)  - 


(3.8) 
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Here  w^1^,  w*2*,  and  w*3* 


denote  the  components  of  w.  These  are  given 

.(2) 


explicitly  in  (2.3).  We  note  that  (3.8)  is  a  perturbation  of  wv 


(|) 

This  expression  can  be  simplified  somewhat.  The  factor  (1  -  $8 ) 


.-1 


can  be  approximated  by  1  +  o(jjr  <J> )  •  When  this  is  multiplied  by  ( i£ )  the 

result  is  (i£)-1  +  o(~$).  For  reasons  stated  earlier,  terms  of  order 

K 

—  $  can  be  omitted  without  affecting  the  order  of  the  error  terms  in  the 

e2 

system  (3.6).  We  can  therefore  replace  (3.8)  with  a  new  fast  quantity 

si(«y+£t)  ^ 2  ff  («}  a_  _  1  —  '  -  -  -(1) 

.  5  ll  U  oil  U 

'(3) 
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We  need  to  find  a  condition  which  suppresses  (3.9)  at  the  boundary 

x  *  0.  If  the  coefficients  are  independent  of  y  and  t,  then  the  bracketed 

quantity  in  the  integrand  is  the  Fourier  transform  of  (3.9),  give  or  take  a 
2 

factor  $  ,  and  we  can  accomplish  what  we  want  by  setting  this  quantity  equal 
to  zero  at  the  boundary.  If  we  do  this,  clear  denominators,  and  then  invert 


the  Fourier  transform,  we  obtain 


0  for  x  -  0 


(4c  >(Y11  “  Y13  +  Y31  "  Y33)w 


In  the  next  section  we  will  say  more  about  the  validity  of  this  boundary 
condition. 

The  condition  (3.10)  is  written  in  terms  of  the  components  of  the 
vector  w  which  appears  in  the  system  (2.2).  He  now  use  the  definition  (2.3) 
of  w  to  write  the  condition  in  terms  of  the  variables  u,  v,  and  p  in  the 
original  system  (1.1).  The  result  is 


It  <tt+P>  +  a  §7  -  Z  Ia(Y12  *  Y32)  +  y°-c)j v 


(3.11) 
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•  (4r,<Y11  ‘Y13  +Y31  -  Y33,(tt“p)  -0 


for  x  -  0 


4.  Discussion  of  the  boundary  condition* 

In  this  section  we  discuss  the  boundary  condition  (3*10)  and  generate 
some  additional  conditions  which  will  be  needed  for  the  test  problem  presented 
in  Section  5* 

If  the  coefficients  in  the  system  depend  on  y  or  t,  the  derivation  of 
(3.10)  from  (3.9)  is  not  valid.  However,  we  can  show  that  (3.10)  is  still 
useful  in  this  case.  Suppose  that  this  condition  holds,  and  write  it  in  the 
simpler  form 


(4.1) 


3w(2)  a  3w(1)  (1)  (3) 

3t  n  3y  +V  +V  “°  ‘ 


if 


If  we  apply  the  operator  having  the  symbol  “  to  (4.1),  the  result  is 
0  -  •  ^"i£"w(2)  +  "(iw  2-  +  Pi)-W(1)  ♦  f3w(3>] 

-  «♦  V2’  ♦  2=  ♦  ^  ♦  -<p)  ♦  <Kp»]* 

. !  ^ ,  !<,;«)  „  ;<2>  t  ^  ;<»>].* 


(4.2) 


+  o(jjr  ♦  2)w  +  (order  (-2))w  . 

According  to  (3.9)  and  the  definitions  of  Fj  and  Fj  implied  by  (4.1),  the 
integral  in  the  last  line  is  the  approximation  to  the  incoming  fast  part  of 
the  solution.  The  entire  last  line  is  equal  to  zero,  so  this  fast  part  must 


o(2—  + 2)w  +  order  (-2)w 


be  equal  to 


at  x  »  0.  The  incoming  fast  part  is  therefore  small  compared  to  w  for 
large  frequencies  and  for  angles  away  from  tangential  incidence. 

It  may  be  worthwhile  to  compare  the  boundary  condition  (3.11),  or  (3.10), 
with  the  condition 

(4.3)  w*2*  ■  u+p  -  0 

which  was  mentioned  early  in  Section  2.  This  condition  was  derived  from  the 

system  (2.2),  in  which  the  coefficient  of  the  x-derivative  is  a  diagonal 

matrix.  The  newer  condition  (3.11)  is  based  on  the  incoming  fast  variable 

(3.8)  which  was  obtained  from  a  more  extensive  uncoupling  of  the  system.  An 

inspection  of  (3.8)  shows  that  this  variable  can  be  considered  a  perturbation 

of  w<2>  ■  u+p,  so  in  some  sense  (3.11)  is  a  refinement  of  (4.3).  One 

obvious  difference  between  the  two  is  the  presence  of  derivatives  in  (3.11). 

This  is  a  result  of  the  need  to  clear  denominators  in  the  Fourier  transform  of 

the  incoming  fast  part.  The  other  difference  is  the  presence  of  terms  which 

do  not  involve  u+p.  The  term  a  |~  *  result  of  uncoupling  the  leading 

symbol,  and  the  other  terms  in  (3.11)  are  the  result  of  reducing  the  coupling 

caused  by  terms  of  order  zero.  The  term  a^(a-c)  corresponds  to  the  part  of 

the  zero-order  coupling  which  resulted  from  the  prior  uncoupling  of  the 

leading  symbol.  If  we  had  not  carried  out  the  lower-order  uncoupling,  then  the 

9  dv 

boundary  condition  would  have  been  (u+p)  +  <*  —  -  0. 

Up  to  now  we  have  discussed  boundary  conditions  only  for  the  incoming 
fast  part  of  the  solution.  If  the  boundary  x  -  0  is  an  inflow  boundary, 
i.e.,  if  o  <  0  in  (2. 2),  then  it  is  also  necessary  to  specify  a  condition 
for  the  slow  part.  One  possibility  is  to  use  the  system  (2.2)  to  obtain  the 
condition 

(4.4)  w^  ■  v  *  given  function,  for  x  *  0 
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Another  possibility  is  to  try  to  base  a  boundary  condition  on  the  more 
extensively  uncoupled  system  (3.6).  We  could  presumably  prescribe  a  value  for 
the  Fourier  transform  of  the  slow  component  in  (3.6),  clear  denominators,  and 
then  apply  an  inverse  transform  to  obtain  an  inhomogeneous  boundary  condition 
analogous  to  (3.11). 

The  first  approach  suggested  here  is  acceptable,  but  the  second  one  is 
not.  The  use  of  the  cutoff  function  $  means  that  we  have  uncoupled  the 
system  only  on  the  wedge  T  in  the  (&>,(;)  space  which  corresponds  to  rapidly 
moving  waves.  This  is  clearly  no  restriction  when  we  are  seeking  boundary 
conditions  which  suppress  the  incoming  fast  part  of  the  solution.  But  in  the 
present  case  it  is  a  major  restriction,  since  the  slow  part  of  the  solution  is 
associated  with  the  entire  (<*>,£)  space.  The  partially  uncoupled  system 
(3.6)  cannot  give  a  full  description  of  the  slow  part,  so  there  is  no  point  in 
trying  to  use  this  system  to  find  an  improvement  of  the  condition  (4.4). 

We  therefore  prescribe  the  conditions 

U  (u+p)  +  -  ta<Y12+Y32)  +  «y(«-c)]v 

14-51  -  lSr)<7ir1,i3-,,'3.-ir33>(“-*'>  * 0 

v  ”  given  function,  for  x  -  0  , 

if  x  -  0  defines  an  inflow  boundary,  in  the  case  of  an  outflow  boundary 
only  the  first  condition  should  be  used. 

Up  to  this  point  the  discussion  has  been  limited  to  problems  defined  on 
the  half-plane  x  >  0.  In  the  next  section  we  present  the  results  of  some 
numerical  computations  for  the  system  (1.1)  defined  on  the  unit  square 
0<x<  1,  0  <  y  <  1.  For  these  computations  we  need  boundary  conditions 
analogous  to  (4.5)  for  the  boundary  segments  y  »  0,  x  *  1,  and  y  =  1. 
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Thai*  conditions  can  bs  obtained  through  suitable  orthogonal  transformations 
of  the  independent  space  variables  x  and  y ,  together  with  the  correspond¬ 
ing  transformations  of  the  velocity  components  u  and  v.  The  necessary 
calculations  are  outlined  in  [1]*  The  results  aret 

Boundary  segment  y  -  Ot 

(y*p>  ♦  »  ♦  #xts-o>w 

-  (5r»" >)><’-»>  *  0 

u  -  given  function. 

Boundary  segment  x  -  1* 

<-u+p)  -  «  ^  "  “yt"0"®*!* 

<«•*>  -  <'^a>«T11-wr1j-Tj1-Tjj)(-»-p)  -  0 

v  •  given  function  . 

Boundary  augment  y  «  1i 

!l  - 1  |a  -  \  t-»(-T21n31)  -  ex(-s-c)i« 

-  <r|s£,<1f2a*T23-T32 -  0 
u  ■  given  function  * 

For  each  segment  the  first  condition  is  the  one  which  is  intended  to 
suppress  the  incoming  fast  part  of  the  solution.  The  second  condition 
prescribes  a  value  for  the  slow  variable,  and  it  should  be  imposed  only  when 
the  segment  is  an  inflow  boundary. 


5.  Numerical  Computations 

In  this  section  we  present  the  results  of  some  numerical  computations 
Involving  the  boundary  conditions  which  have  just  been  derived.  He  consider 
the  system 


on  the  unit  square  0<x<1,  0<y<1.  This  system  Is  a  special  case  of  the 
system  (1.1).  Two  different  choices  for  the  matrix  (y^)  will  be 
considered. 

We  wish  to  compare  three  types  of  boundary  conditions  for  this  system. 

The  first  of  these  is  obtained  by  diagonalizing  the  coefficient  of  the  normal 
derivative  and  then  defining  boundary  conditions  in  terms  of  the  dependent 
variables  in  the  new  system.  These  variables  will  be  referred  to  as 
"characteristic  variables".  This  was  discussed  early  in  Section  2, 
immediately  after  equation  (2.3).  For  the  four  sides  of  the  unit  square  the 
incoming  fast  characteristic  variables  are  the  quantities  in  (4.5)  and  (4.6) 
which  are  differentiated  with  respect  to  time.  The  second  set  of  conditions 
is  obtained  by  uncoupling  the  leading  symbol  in  the  manner  described  earlier, 
but  then  doing  nothing  about  the  zero-order  coupling  in  the  system.  These 
conditions  can  be  obtained  by  deleting  the  zero-order  terms  in  the  derivative 
conditions  appearing  in  (4.5)  and  (4.6).  The  third  set  of  boundary  conditions 
is  obtained  by  also  uncoupling  the  zero-order  terms  in  the  system  of 
differential  equations.  These  are  the  conditions  (4.5)  and  (4.6). 

We  present  two  separate  tests  of  these  conditions,  one  to  demonstrate  the 
effect  of  uncoupling  the  leading  symbol,  and  the  other  to  demonstrate  the 
effect  of  uncoupling  both  the  leading  symbol  and  the  zero-order  term.  In  the 


first  case  we  let  Y^  ■  0  for  all  i,  j  and  use  the  first  two  sets  of 
boundary  conditions*  In  the  second  case  we  use  all  three  sets  of  conditions, 
and  we  let  (Y^ )  be  the  matrix 

(0  10  0 

•10  0  0 

0  0  0 

In  the  computations  we  set  the  solution  equal  to  zero  when  t  -  0.  At 
the  boundary  x  ■  0  we  set  v  (see  (4.5))  equal  to  a  pulse  consisting  of 
half  a  sine  wave  in  t  multiplied  by  half  a  sine  wave  in  the  tangential 
variable  y.  We  use  homogeneous  conditions  on  the  other  boundaries.  The 
nonsero  part  of  the  solution  is  due  entirely  to  the  nonzero  data  at  the 
boundary  x  «  0,  so  it  is  possible  to  study  the  influence  of  these  data  by 
examining  the  size  of  the  solution  in  various  parts  of  the  (x,y)  plane  at 
various  times. 

In  the  computation  the  system  is  approximated  by  the  leap  frog  difference 
scheme.  The  derivative  boundary  conditions  in  (4.5)  and  (4.6)  are 
approximated  by  centered  differences  in  the  time  and  tangent  variables.  The 
outgoing  fast  characteristic  variables  are  extrapolated  at  the  boundary  using 
the  given  differential  equation.  For  this  we  use  centered  differences  in  the 
time  and  tangent  variables,  and  we  approximate  the  normal  derivative  with  a 
forward  difference  which  uses  a  time  average  at  the  back  point.  At  an  outflow 
boundary  the  slow  characteristic  variable  is  extrapolated  in  the  same  manner. 

The  boundaries  y  -  0  and  y  -  1  are  characteristic  for  the  system 
(5.1).  At  these  boundaries  we  integrate  the  slow  characteristic  variable  in 
the  boundary  using  a  centered  difference  approximation.  This  is  an  experiment 
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to  see  if  the  incoming  fast  modes  can  be  activated  at  a  characteristic 
boundary,  in  our  earlier  discussion  we  always  assumed  that  the  boundary  was 
noncharacteristic • 

The  surface  pictured  in  Figure  5. 2  and  5.3  are  graphs  of 
(u2  +  v2  +  p2)^2  as  functions  of  x  and  y  for  fixed  t.  The 
configuration  is  shown  in  Figure  5.1. 


We  show  solutions  at  times  t  *  .125,  .25,  and  .375.  The  fast  mode  entering 
through  the  boundary  x  *  0  has  normal  velocity  4  since  o  *  -1  and 
c  =*  -3.  Pulses  entering  on  this  mode  should  therefore  be  visible  near  the 
nearest  boundary  (x  *  1)  in  the  graphs  for  t  ”  .25. 

In  Figure  5.2  we  show  the  effect  of  uncoupling  the  leading  symbol.  In 
this  case  =  0  for  all  i,j.  Figure  5.2(a)  shows  the  solution 

corresponding  to  the  boundary  conditions  defined  in  terms  of  characteristic 
variables.  The  solution  in  Figure  5.2(b)  corresponds  to  the  more  refined 
boundary  conditions.  The  second  set  of  conditions  is  clearly  more  effective 
at  suppressing  the  incoming  fast  part  of  the  solution. 
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In  Figure  5*3  we  show  tha  affact  of  uncoupling  both  the  Loading  symbol 
and  tha  term  of  order  zero.  In  this  case  the  matrix  (Y^)  ia  given  by 
(5.2).  Tha  simplest  boundary  conditions  are  uaad  in  part  (a).  In  part  (b) 
use  tha  boundary  conditions  obtained  by  uncoupling  tha  leading  symbol  only. 
Tha  boundary  conditions  for  part  (c)  are  obtained  by  uncoupling  both  tha 
leading  symbol  and  tha  term  of  order  saro.  Tha  third  sat  of  conditions  is 
clearly  tha  most  affective. 


Figure  5.2.  Effect  of  uncoupling  the  leading  symbol 


Figure  5.3.  Effect  of  uncoupling  both  the  leading  symbol 
and  the  term  of  order  zero. 

(a)  Solution  using  boundary  conditions  based  A  j  I  |  i 
on  characteristic  variables.  /I  l  I  l  }J- 

JJ  la, 


t  ■  .125 


t  «  .25 


r~- 4  ,'nl 


t  ■  .375 


(b)  Solution  using  boundary  conditions  based 
on  uncoupling  the  leading  symbol. 


t  -  .375 


APPENDIX 


PROPERTIES  OF  PSEUDO-DIFFERENTIAL  OPERATORS 
In  this  appendix  we  define  pseudo-differential  operators  and  state 
without  proof  of  some  of  their  basic  properties.  More  extensive  treatments 
can  be  found  in  Nirenberg  [3] ,  Taylor  [4] ,  and  Treves  [5] . 

We  must  first  establish  some  notation.  Partial  derivatives  in  tf1  will 

ct 

be  denoted  by  «  (  where  a  =  (a  , . . . ,a  ) ,  and 

1  n 

„  a.  o  *  a.  a  <* 

3  -31  —  3n“-  <3T>  n  • 

1  n 

The  are  nonnegative  integers.  Differential  operators  can  then  be  written 


in  the  form 


P  -  I  %3‘ 


The  coefficients  aQ  are  functions  on  R°,  and  the  sum  is  taken  over 

finitely  many  multi-indices  a.  we  will  allow  the  possibility  that  P  may 

act  on  vector-valued  functions.  In  that  case  the  a_  may  be  either  scalars 

a 

or  matrices. 


The  Fourier  transform  on  it*1  will  be  denoted  by 
u(£)  -  w  e"ix’5u(x)dx  , 


(  2K  )"  Rn 


where  x  •  K  *  xj^j*  The  ^nverse  Courier  transform  is  then  given  by 

/  eix’5  u(5)d5  . 

r" 

Differential  operators  can  be  represented  in  terms  of  the  Fourier 
transform.  For  suitable  functions  u,  we  have 


(Pu)(x)  a  Z  aa(x)3  u(x) 


Z  a  (x>3 ®  /  eix’5  u 


(5)d€ 


/  eix‘£ 


[Z  a fl|(x)(U)  ]u(C)dC  . 


“l  °n 

Here  (15)°  denotes  the  product  { i€ 1 )  ’  (iCn>  "•  The  equation  can  be 

written  in  the  form 


(A.1) 


(Pu)(x)  -  /  eix  ^  p(x,C)u(C)dC  , 


where  p(x,C)  -  £  aa(x)(i£)  •  The  function  p  is  sometimes  called  the  symbol 
of  the  operator  P. 

Pseudo-differential  operators  are  obtained  by  allowing  a  larger  class  of 
symbols  to  be  used  in  (A. 1).  Every  differential  operator  is  a  pseudo-differential 
operator,  but  not  vice  versa.  One  fairly  general  symbol  class  is  the  class 

Sm  ,,  0  <  6  <  p  <  1,  which  was  introduced  by  HSrmander.  This  is  defined  to 
P#« 

be  the  set  of  all  C  functions  p  which  satisfy  estimates  of  the  form 
(A. 2)  |3®3^p(x,5)|  <  <*  ^(1  +  |Cj)»-Pl®l+6iei  ,  x  €  K,  C  €  Rn 

for  all  a,  3  and  for  every  compact  subset  K  of  R11 .  The  constant  is 
allowed  to  depend  on  a#  3,  and  K.  The  symbol  of  a  differential  operator  of 
order  m  having  smooth  coefficients  clearly  belongs  to  the  class  S™  For 

the  operators  considered  in  this  paper  we  always  have  p  *  1  and  6  ■  o.  In 

general,  the  number  m  appearing  in  (A. 2)  is  called  the  order  of  the 

operator  P  whose  symbol  is  p.  The  order  need  not  be  positive,  and  it  need 

not  be  an  integer. 

If  u  €  CQ,  then  Pu  €  C  .  It  is  possible  to  extend  P  so  that  Pu  is 
defined  for  any  distribution  u  having  compact  support.  In  this  case  Pu  is 
a  distribution. 

If  an  operator  P  has  negative  order,  then  Pu  is  smoother  than  u, 

A 

since  u(C)  is  multiplied  by  a  symbol  which  tends  to  zero  as  1C  I  ♦  ".  P  is 
said  to  have  order  -  •  if  P  €  Sra  for  all  ra  <  0. 
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A  useful  concept  Is  that  of  an  asymptotic  expansion  of  a  symbol*  Suppose 


that  (mj}  j.Q  is  *  sequence  of  real  numbers  such  that  otj  >  for  all 

j  and  m .  ♦  -  •  as  j  ♦  •*  Let  {p.}  be  a  sequence  of  symbols  such  that 
m  3  3 

p^  €  Sp Jg  for  each  j.  A  symbol  p  Is  said  to  be  an  asymptotic  sum  of  the 
Pj ,  written 


P  ~  I  P4  ' 

j-0  3 

k  ®v+i 

provided  p  -  e  sp  6  toT  a3,3,  11  *  T^at  1  s,  the  error  In  each  partial 

sum  must  have  the  same  order  as  the  first  term  omitted  from  the  partial  sum. 
This  concept  Is  analogous  to  the  usual  concept  of  asymptotic  expansion.  In 
fact,  If  a  function  p{£ )  of  one  variable  has  an  asymptotic  expansion 

p(C)  -  l  4  as  {  ♦  • 
j— m  5  3 

In  the  usual  sense,  then  this  expansion  is  also  asymptotic  in  the  sense 
described  above. 

Pseudo-differential  operators  can  be  composed.  Let  P  and  Q  be 
operators  with  symbols  p(x,€)  and  q(x,€),  respectively,  and  suppose  that 
q  has  compact  support  in  x.  The  composition  P(Qu)  is  then  well-defined 
and  is  given  by  a  pseudo  differential  operator  whose  symbol  has  the  asymptotic 
expansion 

(A. 3)  0  ~  l  -r-j —  3jp(x,5)»“q(x,?)  . 

*0  |o  |>0  i,a|ai  4  x 

The  sum  is  taken  over  all  multi-indices  a  ”  (a^,...,a  )  having  nonnegative 

components.  The  order  of  a  is  given  by  |a|  *  Ea^,  and  the  factorial  a! 

denotes  the  product  a „  la  *...*  a  I. 

12  n 

It  follows  from  (A. 2)  that  when  a  symbol  is  differentiated  with  respect 
to  C,  the  result  is  a  symbol  of  lower  order.  This  implies  that  the  leading 
order  term  in  (A. 3)  corresponds  to  |a|  *  0  and  is  equal  to  p( x,£ )q( x,5 ) . 
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The  symbol  of  the  product  of  two  operators  is  therefore  equal  to  the  product 

of  their  symbols,  up  to  certain  terns  of  lower  order. 

The  makes  sense  when  we  consider  the  special  case  of  differential 

.  a  ,  8 

operators.  The  composition  of  two  operators  a(x)3  and  b(x)d  is  equal  to 
,a  .  8  a+B 

ad  (bd  )  “  abd  +  lower  order  terms  involving  derivatives  of  b.  In  this 

case  the  composition  law  can  be  derived  using  Leibniz'  rule.  For  general 

pseudo-differential  operators  the  derivation  is  much  more  complicated. 

In  this  paper  we  use  the  concept  of  parametrix,  or  approximate  inverse, 

of  a  pseudo-differential  operator.  A  parametrix  of  an  operator  Q  is  an 

operator  such  that  QQ~ 1  -  I  and  Q~^Q  -  I  are  operators  of  order 

-  **.  At  high  frequencies  Q~ 1  is  very  nearly  a  true  inverse  of  Q.  Q~ 1  is 

uniquely  determined  up  to  order  -  <*. 

Let  q  be  the  symbol  of  the  operator  Q.  If  the  function  q  is  a 

nonzero  scalar  or  a  nonsingular  matrix,  then  it  can  be  shown  that  Q-1  exists 

and  that  its  leading  symbol  is  q- 1 .  The  proof  consists  of  writing  a  formal 

expansion  for  Q~ 1  and  solving  for  the  terms  one  by  one. 

An  operator  of  the  form  I  +  K,  where  K  has  negative  order,  has  a 

2 

parametrix  which  is  given  by  the  expansion  (I+K)  -  I  -  K+K  -  ...  . 


I 


I 
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